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We consider a self-consistent system of Bianchi type-I (BI) gravitational field and 
a binary mixture of perfect fluid and dark energy given by a cosmological constant. 
The perfect fluid is chosen to be the one obeying either the usual equation of state, 
i.e., p = £e, with £ G [0, 1] or a van der Waals equation of state. Role of the A term 
, in the evolution of the BI Universe has been studied. 
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In view of its importance in explaining the observational cosmology many authors have 
considered cosmological models with dark energy. In a recent paper Kremer [lj has modelled 
the Universe as a binary mixture whose constitutes are described by a van der Waals fluid and 
by a dark energy density. Zlatev et al. |2j showed that "tracker field" , a form of qiuntessence, 
may explain the coincidence, adding new motivation for the quintessence scenario. The 
fate of density perturbations in a Universe dominated by the Chaplygin gas, which exhibit 
negative pressure was studied by Fabris et al. [3] . Model with Chaplygin gas was also studied 
in the Refs. 0, H|. In doing so the author considered a spatially flat, homogeneous and 
isotropic Universe described by a Friedmann- Robert son- Walker (FRW) metric. Since the 
theoretical arguments and recent experimental data support the existence of an anisotropic 
phase that approaches an isotropic one, it makes sense to consider the models of Universe 
with anisotropic back-ground in presence of dark energy. The simplest of anisotropic models, 
which nevertheless rather completely describe the anisotropic effects, are Bianchi type-I (BI) 
homogeneous models whose spatial sections are flat but the expansion or contraction rate is 
direction-dependent. In a number of papers, e.g., jH0|' we nave studied the role of a A term 
in the evolution of a BI space-time in presence of spinor and/or scalar field with a perfect 
fluid satisfying the equation of state p = (e. In this paper we study the evolution of an 
initially anisotropic Universe given by a BI spacetime in presence of a perfect fluid obeying 
not only p = (e, but also the Van der Waals equation of state. 
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II. BASIC EQUATIONS 

The Einstein field equation on account of the cosmological constant we write in the form 

r;~Is;r = kt; + 8;a. (2.1) 

Here R 1 ^ is the Ricci tensor, R is the Ricci scalar and k is the Einstein gravitational con- 
stant. As was mentioned earlier, A is the cosmological constant. To allow a steady state 
cosmological solution to the gravitational field equations Einstein |s| introduced a funda- 
mental constant, known as cosmological constant or A term, into the system. Soon after E. 
Hubble had experimentally established that the Universe is expanding, Einstein returned to 
the original form of his equations citing his temporary modification of them as the biggest 
blunder of his life. A term made a temporary comeback in the late 60's. Finally after the 
pioneer paper by A. Guth [^| on inflationary cosmology researchers began to study the mod- 
els with A term with growing interest. Note that in our previous papers |y, |3] we studied 
the Einstein field equations where the cosmological term appears with a negative sign. Here 
following the original paper by Einstein and one by Sahni 0| we choose the sign to be 
positive. In this paper a positive A corresponds to the universal repulsive force, while a 
negative one gives an additional gravitational force. Note that a positive A is often taken 
to a form of dark energy. 

We study the gravitational field given by an anisotropic Bianchi type I (BI) cosmological 
model and choose it in the form: 

ds 2 = dt 2 - a 2 dx 2 - b 2 dy 2 - c 2 dz 2 , (2.2) 

with the metric functions a, b, c being the functions of time t only. 

The Einstein field equations (j2.1J) for the BI space-time in presence of the A term now 
we write in the form 
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Here over-dot means differentiation with respect to t. The energy-momentum tensor of the 
source is given by 

t; = (e+ p )uX-ps;, (2.4) 

where VtF" is the flow vector satisfying 

g^u» = l. (2.5) 

Here e is the total energy density of a perfect fluid and/or dark energy density, while p is the 
corresponding pressure, p and e are related by an equation of state which will be studied 
below in detail. In a co-moving system of coordinates from ()2.4|) one finds 

T ° = e, Tl = T 2 = Ti = -p. (2.6) 
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In view of ()2.6|) from ()2.3|) one immediately obtains 

/dt 

b(t) = D 2 r^ 3 e W [X 2 J -^-], (2.7b) 

c(t) = ZV 1/3 exp[X 3 J A]. (2.7c) 

Here Di and Xi are some arbitrary constants obeying 

D X D 2 D 3 = 1, X 1 +X 2 + X 3 = 0, 
and t is a function of t defined to be 

t = abc. (2.8) 

From (J2.3j) for r one find 

T - = ^(e- P )+3A. (2.9) 
On the other hand the conservation law for the energy-momentum tensor gives 



i= — (e + p). (2.10) 

T 

After a little manipulations from ()2.9|) and ()2.10|) we find 

f 2 = 3(ke + A)t 2 + d, (2.11) 

with Ci being an arbitrary constant. Let us now, in analogy with Hubble constant, define 

t a b c , 
- = - + - + - = 3H. 2.12 

t a b c 

On account of ()2.12j) from (|2.11[) one derives 

ke = 3H 2 - A-d/(3r 2 ). (2.13) 

It should be noted that the energy density of the Universe is a positive quantity. It is 
believed that at the early stage of evolution when the volume scale r was close to zero, the 
energy density of the Universe was infinitely large. On the other hand with the expansion 
of the Universe, i.e., with the increase of r, the energy density e decreases and an infinitely 
large r corresponds to ae close to zero. Say at some stage of evolution e is too small to be 
ignored. In that case from ()2.13j) follows 

3# 2 -A^0. (2.14) 

As it is seen from (I2.14J) in this case A is essentially non-negative. We can also conclude 
from 1)2.14)1 that in absence of a A term beginning from some value of r the evolution of 
the Universe comes stand-still, i.e., r becomes constant, since H becomes trivial, whereas in 
case of a positive A the process of evolution of the Universe never comes to a halt. Moreover 
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it is believed that the presence of the dark energy (which can be explained with a positive 
A as well) results in the accelerated expansion of the Universe. As far as negative A is 
concerned, its presence imposes some restriction on e, namely, e can never be small enough 
to be ignored. It means in that case there exists some upper limit for r as well (note that 
t is essentially nonnegative, i.e. bound from below). In our previous papers we came to the 
same conclusion @,@] [with a positive A which in the present paper appears to be negative]. 
Inserting (f!H!?J) and (EH3|) into $Bj\) one now finds 



H = _l (3H «_A + ^ + ^) = -5( e+P )-^. (2.15) 

In view of (12. 13(1 from (j2.15j) follows that if the perfect fluid is given by a stiff matter where 
p = e, the corresponding solution does not depend on the constant C\. 

Let us now go back to the Eq. (J2.11j) . It is in fact the first integral of (J2.9|) and can be 
written as 

f = ±y/C x + 3(/«£ + A)r 2 (2.16) 
On the other hand, rewriting (j2.10|) in the form 



T 



e + p t 1 



(2.17) 



and taking into account that p is a function of e, one concludes that the right hand side of 
the Eq. ()2.9|) is a function of r only, i.e., 

f = y(e-j9)r + 3Ar = ^(r). (2.18) 

From a mechanical point of view Eq. (|2.18|) can be interpreted as an equation of motion 
of a single particle with unit mass under the force F{t). Then the following first integral 
exists 0: 

t=y/2[£-U(r)]. (2.19) 

Here £ can be viewed as energy and U{t) is the potential of the force T . Comparing the 
Eqs. flUEl and flOU) one finds 8 = C x /2 and 

U(r) = -^(ke + A)t 2 . (2.20) 

Let us finally write the solution to the Eq. (|2.9|) in quadrature: 

dr 



y/Ct + 3(ks + A)r 2 



t + t , (2.21) 



where the integration constant t can be taken to be zero, since it only gives a shift in time. 

In what follows we study the Eqs. ()2.9j) and (|2.1()j) for perfect fluid obeying different 
equations of state. 



III. UNIVERSE FILLED WITH PERFECT FLUID 

In this section we consider the case when the source field is given by a perfect fluid. Here 
we study two possibilities: (i) the energy density and the pressure of the perfect fluid are 
connected by a linear equation of state; (ii) the equation of state is a nonlinear (Van der 
Waals) one. 
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A. Universe as a perfect fluid with p p f = (e p { 

In this subsection we consider the case when the source field is given by a perfect fluid 
fluid obeying the equation of state 

P P f = C^pf, (3-1) 

where ( is a constant and lies in the interval C G [0, 1]. Depending on its numerical value, 
£ describes the following types of Universes 

( = 0, (dust Universe), (3.2a) 

( = 1/3, (radiation Universe), (3.2b) 

C G (1/3, 1), (hard Universes), (3.2c) 

( = 1, (Zel'dovich Universe or stiff matter). (3. 2d) 

In view of (|3.1j) . from (j2.1()|) for the energy density and pressure one obtains 

e pi = e /r^\ p pi = Ce /r^\ (3.3) 
where Eq is the constant of integration. For r from (|2.21j) one finds 

dr 



y/d + 3(^ tW + At 2 ) 



t. (3.4) 



As one sees, the positivity of the radical in ()3.4|) for a negative A imposes some restriction 
on the upper value of r, i.e., r should be bound from above as well. In Fig. ^the graphical 
view of the potential U(t) is illustrated for a negative A. As it was mentioned earlier, £ or 
C\ in case of ( = 1 does not play any role. Universe in this case initially expands, reaches 
to the maximum and then begin to contract finally giving rise to a space-time singularity [cf 
Fig. E]. For the other cases depending on the choice of £ expansion of the Universe is either 
non-periodic [Fig. Ej with a singularity at the end or oscillatory one without space-time 
singularity [Fig. In Fig. 0] we demonstrate the evolution of the BI Universe with a 
positive A. In this case the Universe expands exponentially, the initial anisotropy quickly 
dies away and the BI Universe evolves into a isotropic FRW one. There does not any upper 
bound for r in case of a positive A. Note that in the Figs. (P-EJ) d, r, h and s stand for 
dust, radiation, hard Universe and stiff matter, respectively. 



In absence of the A term one immediately finds 

T = AtW+V, (3.5) 

with A being some integration constant. As one sees from ()2.7|) . in absence of a A term, 
for C < 1 the initially anisotropic Universe eventually evolves into an isotropic FRW one, 
whereas, for ( = 1, i.e., in case of stiff matter the isotropization does not take place. 



B. Universe as a van der Waals fluid 



Here we consider the case when the source field is given by a perfect fluid with a van der 
Waals equation of state in absence of dissipative process. The pressure of the van der Waals 
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FIG. 1: View of the potential W(r). As one 
sees in case of stiff matter this potential al- 
lows only non-periodic solution. 



FIG. 2: Evolution of volume scale r with a 
negative A and C = —0.1. As one sees, in 
this case the model with perfect fluid given 
by dust, radiation and hard Universe allow 
oscillation, whereas, stiff matter gives rise 
to a non-periodic solution. 
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FIG. 3: Evolution of volume scale r with a 
negative A and C = 0. In this case the model 
with perfect fluid given by dust, radiation, 
hard Universe and a stiff matter gives rise to 
a non-periodic solution. 
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FIG. 4: Evolution of the Universe with a 
positive positive A. In this case indepen- 
dent to the choice of zeta the expansion of 
the Universe is always exponential. 



fluid p w is related to its energy density e w by [1 

8We w 2 

Pw = — 3e w . (3.6) 

In (|3.6|) the pressure and the energy density is written in terms of dimensionless reduced 
variables and W is a parameter connected with a reduced temperature. In the Figs. Eland 
E] the energy density and the pressure of the system are illustrated with a negative and a 
positive A term as well as in absence of it. 
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FIG. 5: View of energy density e and pres- 
sure p in case of a Van der Waals fluid with 
a negative A. 
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FIG. 6: View of energy density e and pres- 
sure p in case of Van der Waals fluid with 
A > 0. 
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FIG. 7: Evolution of r(t) with the BI Uni- 
verse filled with Van der Waals fluid. Inde- 
pendent to the sign of A the model provides 
provides with exponentially expanding Uni- 
verse. 



FIG. 8: Evolution of the BI Universe filled 
with different types of perfect fluid in ab- 
sence of a A term. As one sees, in case of 
Van der Waals fluid r(t) grows faster at the 
early stage, then slows down with time. 



Inserting ()3.6|) into (|2.15J) on account of (j2.13j) one finds 



H 



{3H 2 - A - C 1 /(3r 2 )}[(3 + SW)k - {3H 2 - A - CV(3t 2 )}] 
2(3/?-{3# 2 -A-C 1 /(3r 2 )}) 

!-({3tf 2 -A-Cy(3r 2 )}) 2 . 



(3.7) 



It can be easily verified that the Eq. ()3.7)) in absence of A term and C\ = and k = 3 
coincides with that given in : 



H 



H 2 | 8WH 2 



(3.8) 
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The solution of the second-order differential equation (|3.7|) for H(t) can be found by 
specifying the initial value for H(t) at t — 0, for a given value of parameter W. Here we 
graphically present some results concerning the evolution of BI Universe with a Van der 
Waals fluid. In Fig. [7| we compare the evolution of r with and without A term. As one 
sees, the character of evolution does not depend on the sign of A. In all cases we find 
exponentially expanding Universe, though the rapidity of growth depends on A. The Fig. |H] 
gives the comparison of the expansion of r with perfect fluid obeying different equations of 
state. 



IV. CONCLUSION 

The evolution of an anisotropic Universe given by a Bianchi type I cosmological model 
is studied in presence of a perfect fluid and a A term. It has been shown that in case of 
a perfect fluid obeying p = (e, where p and e are the pressure and energy density of the 
fluid, respectively, a negative A may generate an oscillation in the system thus giving rise 
to a singularity-free mode of expansion. Introduction of a positive A in this case results 
in a rapid expansion of the Universe. If the Universe is filled with a Van der Waals fluid, 
no oscillatory or non-periodic mode of expansion occurs. Independent to the sign of A the 
Universe in this case expands exponentially. 
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